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Recent observations of type la supernovae indicate that the Universe is in an accelerating phase 
of expansion. The fundamental quest in theoretical cosmology is to identify the origin of this 
phenomenon. In principle there are two possibilities: 1) the presence of matter which violates 
the strong energy condition (a substantial form of dark energy), 2) modified Friedmann equations 
(Cardassian models - a non-substantial form of dark matter). We classify all these models in terms 
of 2-dimensional dynamical systems of the Newtonian type. We search for generic properties of the 
models. It is achieved with the help of Peixoto's theorem for dynamical system on the Poincare 
sphere. We find that the notion of structural stability can be useful to distinguish the generic cases 
of evolutional paths with acceleration. We find that, while the ACDM models and phantom models 
are typical accelerating models, the cosmological models with bouncing phase are non-generic in the 
space of all planar dynamical systems. We derive the universal shape of potential function which 
gives rise to presently accelerating models. Our results show explicitly the advantages of using a 
potential function (instead of the equation of state) to probe the origin of the present acceleration. 
We argue that simplicity and genericity are the best guide in understanding our Universe and its 
acceleration. 
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As indicated by observations of distant type la supernovae [J, Q and other complementary observations 0, Q our 
Universe is presently accelerating. The problem of accelerated expansion of the current universe seems to be one of 
the most fundamental problems of theoretical physics of the XXI century (for a recent review see 0, @ and references 
, therein). 

' There are several propositions of explanation of the observational fact that the Universe is speeding up rather 
than slowing down. They involve the cosmological constant, time dependent vacuum energy, dynamical scalar field 
(quintessence) or modified Friedmann equations 0,00- However, in all truth, we do not know what is causing 
the effect. While the proposition of the cosmological constant is an attractive idea from both the theoretical and 
observational point of view, it entails a crucial problem. Why is its value, as measured by distant SNIa observations, 
so small when compared to the value obtained from the quantum field theory considerations? Neither do we understand 
^ , why the dark energy density is of the same order of magnitude as the matter density during the present epoch. 

Recently some new developments in dark energy studies were proposed. For example non-Abelian Einstein-Born- 
rS Infeld dilaton theory was explored [H EH, as well as cosmology based on its generalization [l^ . 

In this work, we start form general relativity cosmological models with Robertson-Walker symmetry rather than 
Euler's and Poison's equations. We adopt the particle-like description of FRW cosmologies which we extend to a 
general class of FRW cosmologies with dark energy component parameterized by the scale factor (or redshift). The 
main advantage of such formulation is that the one dimensional potential function contains the complete information 
about the dynamics. Finally, we obtain a unified description of a large class of cosmological models in the notion of a 
particle moving in a one-dimensional potential well. Hence, the evolution of the system can be naturally reduced to 
a simple 2-dimensional dynamical system of the Newtonian type. 

The problem of structural stability was studied by Lazkoz in the context of cosmic acceleration • The ideas 

of rigidity and fragility of solutions are translated from the language of structural stability which was introduced by 
Andronov and Pontryagin |lq . Following Lidsey |l7j the concepts of rigidity and fragility should be described through 
a condition on the functional form of the Hubble function. In this paper we formulate the condition for structural 
stability in terms of the potential function rather than the Hubble function. However, the simple relation between 
them makes the equivalence of this approaches. 

All attempts at cosmological modelling involve at the very beginning a large number of simplification and theoretical 



•Electronic address: uoszydlo@cyf-kr.edu.pl 



2 



assumption regarding the parameters of the dark energy of unknown form. The question then arises whether such 
simple models are representing properties of real systems. Usually, for dynamical systems to be viable as models 
they need to be structurally stable, i.e. their dynamics must preserve its qualitative characteristics under small 
perturbation |l8lll9| |. Such a point of view resolves the "approximation problem" because our models of reality are, 
by definition, not precise, but might be very close to it. Moreover, the structural stability framework can be useful 
from the methodological point of view because a'priori we do not know the functional forms of dynamical systems. 
Of course there is no objective reason to believe that all dynamical systems should be structurally stable. From the 
philosophical view of point, it is a presumption of convenience (or economy). According to Occam's Razor, simple 
theories are more economical and mostly cases a little improvement of prediction is paid a very large increase in 
complexity of a model. To discriminate between different dark energy models, we postulate a simplicity principle that 
among all dynamical laws describing the cosmological evolution, the laws with the smallest complexity are chosen. 
It is interesting that Occam's Razor became a cornerstone of modern theory of induction [20I Ell ]. The notion of 
structural stability was discussed in the context of Kaluza-Klein theories [22I . In this class of models it is required 
the existence of the mechanism of dynamical reduction of extra-dimensions, i.e. the configuration of FRW x {static 
internal space} should be an attractor. 

In this paper we apply the structural stability notion as the discriminatory among the cosmological models with 
acceleration phase of expansion. In other words structural stability acting as Occam's Razor allows us to choose the 
simplest (the ACDM or phantom CDM) cosmological models with acceleration. From a physical view of point these 
two are two-phased models with decelerating matter-dominated phase and accelerating dark energy dominated phase. 

It is interesting that the strong energy condition (p + 3p) > and (p + p) > (SEC) which determines qualitative 
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dynamics of the FRW cosmological models can be translated into the constraint on the potential function 

Ma) = 4t-^- > and 7 v(«) < 2 
dm a 

where Iy (a) is an elasticity coefficient which measures the logarithmic slope of the potential function with respect to 
the scale factor a. 

The famous Hawking- Penrose singularity theorems invoke this condition |2^|. Note that the null energy condition 
(NEC) p + p > involving a derivative of the potential function can be formulated in an exact form 

-^ + 2-^->0, Via) < Aa 2 , A = const > 0. 
da a 

The violation of the SEC means that V(a) is a decreasing function of the scale factor. The violation of the NEC 
denotes that in some region the diagram of the potential function is over the diagram of a parabola V(a) oc a 2 . 

The idea of obtaining a Newtonian analogy to the FRW cosmology within the framework of classical mechanics 
has been considered since the classical papers of Milne and McCrea |2(| [27I • In their formulation of the cosmological 
problem the fluid dynamics for a dust filled universe is applied but there arises a crucial difficulty because pressure p 
does not play a dynamical role like in general relativity theory. As a consequence, there is no classical analogy of a 
radiation filled universe. 

If we assume the validity of the Robertson- Walker symmetry for our universe which is filled with perfect fluid 
satisfying the general form of the equation of state p e ff = w e ff(a)p e ff, then p c g = p c g(a) i.e. both the effective energy 
density p c s and pressure p D f{ are parameterized by the scale factor as a consequence of the conservation condition 

Pcff = -3ff (p off + p cS ), (1) 

where dot denotes differentiation with respect to the cosmological time t and H — (In a) is the Hubble function. 

To unify two main approaches to explaining SNIa data, i.e. 1) dark energy, and 2) modification of FRW equations, 
we generalize the acceleration equation 

- = -]:{p e s + 3p eS ) + -a m , (2) 
a b b 

where B and m are constants, rather than doing so for the first Friedmann equation like in Freese and Lewis's 
approach. 

We assume that the Universe is filled with standard dust matter (together with dark matter) and dark energy X 

Pctt = + W X PX, /ox 

Pcff = Pm + Px, 

where wx = wxifl) is the coefficient of the equation of state for dark energy parameterized by the scale factor or 
redshift z: 1 + z = a" 1 . 
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One can check that the Raychaudhuri equation (JSJ can be rewritten to the form analogous to the Newtonian 
equation 



dV 

da ' 



(4) 



if we choose the following form of the potential function V{a) 



~ l 7 {p c xa 2 + — ^-a m+2 ) for m + -2, 
6 m + 2 

— g (Pcfffl 2 + B In a) for m = — 2, 



(5) 



where p e ff satisfies the conservation condition As the alternative method to obtain JSJ is integration by parts eq. 
Q) with the help of conservation condition QJ. 

If we put B = in formula (JSJ) then we obtain the standard cosmology with dark energy. If we consider m 7^ —2, 
Pcff = Pm, then the Cardassian cosmology can be recovered provided that p — p(a) is assumed. The case of B ^ and 
m = — 2 represents a new exceptional case which appears as a consequence of the generalized Raychaudhuri equation 
instead of the Friedmann first integral which assumes the following form 

a 2 k B _ 
Pcff-3-=3- — a m , (6) 

i & rt ^ im —I— / 



or 



where 



= -2V (7) 



6 2 6(m + 2) 

Equation Q is the form of the first integral of the Einstein equation with Robertson- Walker symmetry called the 
Friedmann energy first integral. Formally, curvature effects as well as the Cardassian term can be incorporated into 
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the effective energy density (p k = -^pcard 

The form of equation suggests a possible interpretation of the evolutional paths of cosmological models as the 
motion of a fictitious particle of unit mass in a one dimensional potential parameterized by the scale factor. Following 
this interpretation the Universe is accelerating in the domain of configuration space {a: a > 0} in which the potential 
is a decreasing function of the scale factor. In the opposite case, if the potential is an increasing function of a, the 
Universe is decelerating. The limit case of zero acceleration corresponds to an extremum of the potential function. 
The energy conditions were also confronted with SNIa observations . It is shown that all energy conditions seem 
to have been violated in a recent past of evolution of the Universe. 

It is useful to represent the evolution of the system in terms of dimensionless density parameters = pi/{"iH^) : 
where Hq is the present value of the Hubble function. For this aims it is sufficient to introduce the dimensionless scale 
factor x = a/ clq which measures the value of a in the units of the present value ao, and parameterize the cosmological 
time following the rule ( h r: dt|ifo| = dr. Hence we obtain a 2-dimensional dynamical system describing the 
evolution of cosmological models 



d.r 
d7 

dy _ dV 
dr dx 



y, (9a) 

(9b) 



and y 2 /2 + V(x) =0, 1 + z = 2T 1 . Where 



V(X) =~ {n cff X 2 + ft C ard,0X m+2 + «M } » 



ttcff = n m , x- 3 + Q x ,ox- 3{1+wx \ 
for dust matter and quintessence matter satisfying the equation of state px = wxpx, wx — const. 
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TABLE I: The potential functions for different dark energy models. 



The model 



The form of the potential function 



Independent parameters 



Einstein-de Sitter model 
fim.o = 1, ^a,o = 
fifc.o = 



V(x) = -^ox- 1 



Ho 



ACDM model 
f2 m o + S7a.o = 1 



V{x) = -U On.oaT 1 + Ov.oa: 2 



(n m ,o, Hq) 



FRW model filled with n 
noninteracting multi-fluids 

p = wp with dust matter 
and curvature 



v(x) = -u Vmfix- 1 + n k , + Efei n il0 z- s(1+, "« ) 



(f2 m ,0j &A,0, &n,0, Hq) 



FRW quintessence model with 
dust and dark matter X 
wx < — 1 phantom models 



V(x) = — ^\ fi m ,oa; 1 + fl k ,o + &x,ox~ 



(Omfi, Qx,0, Ho) 



FRW model with generalized 
Chaplygin gas [29. 301 



V(l) = -||fim,oa; 1 + n fc , + » C hapLo(As + ^(l+g) ) 1 +° 



(f2 m ,0, ^Chapl.O, Ho) 



FRW models with dynamical 
equation of state for dark energy 
px = wx(a)px and dust 



V(x) = - I jan.oaT 1 + n K o + QxfiX- 1 exp [/* ^^da] J 



(fim.o, f2jc,0j Ho) 



FRW models with dynamical 
equation of state for dark energy 
coefficient equation of state 
Wx — wo + WiZ 



v(z) = -| {n„ >0 (i + «) + njc, (i + z) 1+3 <™°- 



(fbc,o,wo, Wi,-Ho) 



The form Q of a dynamical system opens the possibility of adopting the dynamical systems methods in investi- 
gations of all possible evolutional scenarios for all possible initial conditions. Theoretical research in this area has 
obviously shifted from finding and analyzing particular cosmological solution to investigating a space of all admis- 
sible solutions and discovering how certain properties (like, for example, acceleration, existence of singularities) are 
"distributed" in this space. 

The system (0 is a Hamiltonian one and adopting the Hamiltonian formalism to the admissible motion analysis 
seems to be natural. The analysis can then be performed in a manner similar to that of classical mechanics. The 
cosmology determines uniquely the form of the potential function V(x), which is the central point of the investigations. 
Different potential functions for different propositions of solving the acceleration problem are presented in Tableland 
Table HU 

The key problem of this paper is to investigate the geometrical and topological properties of the multiverse of 
models of accelerating universes which we define in the following way 

Definition 1 By multiverse of models of accelerating Universes we understand the space of all 2- dimensional sys- 
tems of the Newtonian type x — y, y = —dV/dx with suitably defined potential function of the scale factor, which 
characterize the physical model of the dark energy or modification of the FRW equation. 

The organization of the text is the following. In section [H] we investigate the property of structural stability of 
different subsets of multiverse. Section inTl is devoted to investigation of the inverse problem of accelerating cosmology 
- reconstruction of potential function from SNIa data and estimation of the value of transition redshift and Hubble 
function at the moment of changing decelerating phase into accelerating one. In section llVl we distinguish generic 
accelerating models of the multiverse with the help of Peixoto's theorem which characterize the structurally table 
systems on the compact two-dimensional phase plane defined on a two-dimensional Poincare sphere. In section IVI we 
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TABLE II: The potential functions for different cosmological models which offer the possibility of explanation of acceleration 
in terms of modified FRW equations. 



The model 


The form of the potential function 


Independent parameters 


Non-flat Cardassian models 

T~l 1 1 CXf\ r"l"\ T f\ 1 1 TV1 O "f~"f" \i 7<\ 

llllcLl Uy LILlbL IIldlLtl | / . O 




(flo, "m,0, "Card.Oj 


Bouncing cosmological models 

(H/Ho) 2 = n m , x~ m - n n , x~ n , 

/£ ^> III lOl 


v(x) = -| {n m , i- m+2 - n n , x~ n + 2 } 


(Ho, fi m ,o, f2n,o, m,n) 


Randall- Sundrum brane 
models with dust on the brane 
and dark radiation [32] 


V ( x ) = "I {^m.oa:" 1 + M xfi x- 6 + n fc , + Oj,oz~ 4 } 


(Ho, f2m,0, ^A,o) 


Cosmology with spin 

and dust 
(MAG cosmology) [33] 




(-ffo, fis.O, fife.o) 


Dvali, Deffayet, Gabadadze 
brane models 
(DDG) [34] 


V(a;) = -§ wftm.oa; -1 + + \/ n rc,o | 


(-ffo, f2 m ,o) 


Sahni, Shtanov 
brane models [3^, 36j 


V(x) = -ijfim.oa;- 1 + fi CT , + 2fi;,o 
±2^/!T2;,() ^ fim.oZ -1 + f2cr,0 + fij,o + fiAi,,o| 


(-ffo, fim,0, ^CT,0, f2i,0, f^Ajj.o) 


FRW cosmological models 
of nonlinear gravity C oc R n 
with matter and radiation [37, 38] 


H*) = ~h {tt^ox 1 + ^z^i n r , oa; 2 }n nonl , s ( - ! 


(Ho, f2m,0, f^nonl.o) 


ADGP model [3£, 4Q] 
screened cosmological 
constant model 


= -I,; 2 {-^ + ^(2 + +4fi m , (x-3 - I)} 


(r Q H ,H ) 



formulate the conclusions and examine the significance of the obtained results for the philosophical discussion over 
McMullin's indifference principle, and the fine-tuning principle in the modern cosmological context. 

II. STRUCTURAL STABILITY ISSUES 

Einstein's field equations constitute, in general, a very complicated system of nonlinear, partial differential equations, 
but what is made use of in cosmology are the solutions with prior symmetry assumptions postulated at the very 
beginning. In this case, the Einstein field equations can be reduced to a system of ordinary differential equation, 
i.e. a dynamical system. Hence, in cosmology the dynamical systems methods can be applied in a natural way. The 
applications of these methods allow to reveal some stability properties of particular solutions, visualized geometrically 
as trajectories in the phase space. Hence, one can see how large the class of the solutions leading to the desired property 
is, by means of attractors and the inset of limit set (an attractor is a limit set with an open inset - all the initial 
conditions that end up in some equilibrium state) . The attractors are the most prominent experimentally, because of 
the probability for an initial state of the experiment to evolve asymptotically to the limit set being proportional to 
the volume of the inset. 

The idea, now called structural stability, emerged in the history of dynamical investigations in the 1930's with the 
writings of Andronov, Leontovich and Pontryagin in Russia (the authors do not use the name structural stability 
but rather the name "roughly systems" ) . This idea is based on the observation that actual state of the system can 
never be specified exactly and application of dynamical systems might be useful anyway if it can describe features of 
the phase portrait that persist when state of the system is allowed to move around (see Ref. [4ll p. 363] for more 
comments) . 

Among all dynamicists there is a shared prejudice that 
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1. there is a class of phase portraits that are far simpler than arbitrary ones which can explain why a considerable 
portion of the mathematical physics has been dominated by the search for the generic properties. The exceptional 
cases should not arise very often in application and they de facto interrupt discussion (classification) 0, p. 
349]; 

2. the physically realistic models of the world should posses some kind of structural stability because to have many 
dramatically different models all agreeing with observations would be fatal for the empirical methods of science 

El El El El E El m m . 

In cosmology a property (for example acceleration) is believed to be physically realistic if it can be attributed to 
generic subset of the models within a space of all admissible solutions or if it possesses a certain stability, i.e. if it 
is shared by a "epsilon perturbed model" . For example G. F. R. Ellis formulates the so called probability principle 
"The universe model should be the one that is a probable model" within the set of all universe models and a stability 
assumption which states that "the universe should be stable to perturbations" [iij . 

The problem is how to define 

1. a space of states and their equivalence, 

2. a perturbation of the system. 

The dynamical system is called structurally stable if all its 5-perturbations (sufficiently small) have an epsilon 
equivalent phase portrait. Therefore for the conception of structural stability we consider a S perturbation of the 
vector field determined by right hand sides of the system which is small (measured by delta) . We also need a conception 
of epsilon equivalence. This has the form of topological equivalence - a homeomorphism of the state space preserving 
the arrow of time on each trajectory. In the definition of structural stability consider only deformation of the "rubber 
sheet" type stretches or slides of the phase space by a small amount measured by epsilon. 

Fig. H illustrates the property of structural stability of a single spiral attractor (focus) and saddle point, and the 
structural instability of center. The addition of a delta perturbation pointing outward (no matter how week) results 
in a point repellor. We call such a system structurally unstable because phase portrait of the center and focus are 
not topologically equivalent (note that all phase curves around the center are closed in contrast to the focus). Hence 
one can claim that the pendulum system (without friction) is structurally unstable. 

The idea of structural stability attempts to define the notion of stability of differential deterministic models of 
physical processes. 

For planar dynamical systems (as is the case for the models under consideration) Peixoto's theorem Eij states that 
structurally stable dynamical systems form open and dense subsets in the space of all dynamical systems defined 
on the compact manifold. This theorem is a basic characterization of the structurally stable dynamical systems on 
the plane which offers the possibility of an exact definition of generic (typical) and non-generic (exceptional) cases 
(properties) employing the notion of structural stability. Unfortunately, there are no counterparts of this theorem in 
more dimensional cases when structurally unstable systems can also form open and dense subsets. For our aims it is 
important that Peixoto's theorem can characterize generic cosmological models in terms of the potential function. 

While there is no counterpart of Peixoto's Theorem in higher dimensions, it is easy to test whether a planar 
polynomial system has a structurally stable global phase portrait. In particular, a vector field on the Poincare sphere 
will be structurally unstable if there is a non-hyperbolic critical point at infinity or if there is a trajectory connecting 
saddles on the equator of the Poincare sphere S 2 . In opposite case if additionally the number of critical points and 
limit cycles is finite, / is structurally stable on S 2 (see [5l p. 322]). Following Peixoto's theorem the structural 
stability is a generic property of the C 1 vector fields on a compact two-dimensional differentiable manifold A4. 

Let us introduce the following definition 

Definition 2 If the set of all vector fields f G C r (M) (r > 1 ) having a certain property contains an open dense 
subset ofC r (M), then the property is called generic. 

From the physical point of view it is interesting to known whether a certain subset V of C r (M) (representing the 
class of cosmological accelerating models in our case) contains a dense subset because it means that this property 
(acceleration) is typical in V (see Fig. P). 

It is not difficult to establish some simple relation between the geometry of the potential function and localization 
of the critical points and its character for the case of dynamical systems of the Newtonian type: 

1. The critical points of the systems under consideration x = y, y — —^r lie always on the x axis, i.e. they 
represent static universes yo — 0, x — xq; 

2. The point (xo,0) is a critical point of the Newtonian system iff it is a critical point of the potential function 
V(x), i.e. V(x) — E (E is the total energy of the system E — y 2 /2 + V(x); E = for the case of flat models 
and E = —k/2 in general); 
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FIG. 1: Illustration of the structural instability of center and (b) structural stability of saddle point. A perturbation of the 
original vector field / means the addition of a small vector g. The effect of perturbation modifies it at every point in the state 
space (see (a)). The effect of adding a global perturbation to the original system disturbs a center point (see (b)) and preserves 
a saddle point as illustrated with figure (c). 



3. If (xq, 0) is a strict local maximum of V{x), it is a saddle type critical point; 

4. If (xo, 0) is a strict local minimum of the analytic function V(x), it is a center; 

5. If (xo, 0) is a horizontal inflection point of the V(x), it is a cusp; 

6. The phase portraits of the Newtonian type systems have reflectional symmetry with respect to the y axis, i.e. 
x -» x, y -> -y. 

All these properties are simple consequences of the Hartman-Grobman theorem which states that near the non- 
degenerate critical points (hyperbolic) the original dynamical system is equivalent to its linear part. Therefore the 
character of a critical point is determined by the eigenvalues of the linearization matrix given by a simple equation 
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FIG. 2: The model of an accelerating universe given in terms of the potential function and its phase portrait. The domain of 
acceleration is represented by shaded area. Note that only in this case with a single maximum of the potential function there 
are two phases of the universe evolution where the acceleration phase follows the deceleration phase. Note also, that the order 
of appearing of acceleration epochs is important and only in this case we obtain a structurally stable phase portrait. Therefore, 
from Peixoto's theorem we obtain the generic phase portrait for the accelerating universe. It is equivalent to the ACDM model 
scenario. 



A 2 + dct A = 0, where 

A = 



1 

o 

d x J (x ,0) 



is the linearization matrix. Hence, in the case of a maximum we obtain a saddle with Ai, A2 real of opposite signs, 
and if the potential function assumes a minimum at the critical point we have a center with Ai^ purely imaginary of 
mutually conjugate. Therefore, among all distinguished cases, only if the potential function admits a local maximum 
at the critical point we have a structurally stable global phase portrait. Because V < and dV/da = + 3p)(a) 
the Universe is decelerating if the strong energy condition is satisfied and accelerating if the strong energy condition 
is violated. Hence, among all simple scenarios, the one in which deceleration is followed by acceleration is the only 
structurally stable one (see Fig. 01 . 

Let us consider two types of scenarios of cosmological models with matter dominated and dark energy dominated 
phases 

1. the ACDM scenario, where the early stage of evolution is dominated by both baryonic and dark matter, and 
late stages are described by the cosmological constant effects. 

2. the bounce instead initial singularity squeezed into a cosmological scenario; one can distinguish cosmological 
models early bouncing phase of evolution (caused by the quantum bounce) |5l) from the classical bouncing 
models at which the expansion phase follows the contraction phase. In this paper by bouncing models we 
understand the models in the former sense (the modern one). 

For the first class of models we obtain a global phase portrait equivalent to that which is demonstrated in Fig. El The 
eigenvalues of this system are real of opposite signs and this point is a saddle. It describes a stationary but unstable 
universe which is quite similar to the static Einstein Universe filled with both dust matter and with cosmological 
constant. There are, in principle, three representative scenarios of evolution. The trajectories moving in the region B 
confined by the separatrices correspond to the closed universes contracting from the unstable de Sitter node towards 
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FIG. 3: The model of an accelerating universe given in terms of the potential function and its phase portrait. The domain 
of acceleration is represented by shaded area. In this case the existence of two maxima induces the appearance of minima 
of the potential function whose presence gives rise to structural instability. Therefore, from Peixoto's theorem we obtain the 
non-generic phase portrait for the universe accelerating in two domains. 
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FIG. 4: The model of an accelerating universe given in terms of the potential function and its phase portrait. The domain of 
acceleration is represented by shaded area. In this case the universe is accelerating for all trajectories. In the phase portrait 
there is no critical point in the finite domain. The potential function is decreasing function of the scale factor. While this system 
is structurally stable there is no matter dominating phase in the model. And from physical view of point is not interesting. 
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FIG. 5: The model of an accelerating universe given in terms of the potential function and its phase portrait. The domain of 
acceleration is represented by shaded area. The phase portrait of extended bouncing models with the cosmological constant. 
While the early evolution is dominated by matter terms, for late times the cosmological term is dominating. There are three 
characteristic types of evolution: I - inflectional, O - oscillating, B - bouncing. The model is structurally unstable because of 
the presence of a non-hyperbolic critical point (center) on the phase portrait. 

the stable de Sitter node. They are sometimes called bouncing models (it is not the sense of bouncing models in the 
sense used by us in this paper - see Tables III and UTI for comparison). All trajectories on the phase plane are divided by 
a parabolic curve (representing the trajectory of the fiat model) into two disjoint classes of closed and open models. 
The trajectory situated in the region / confined by the upper branch of the trajectory k — and by the separatrices 
(coming out and approaching the saddle) of the saddle point, correspond to the closed universes expanding from 
the initial singularity (x = 0, x = oo) to the stable deS + attractor. Quite similarly, the trajectories located in the 
symmetric region x — > x and y — > — y correspond to the closed contraction universes (y < 0) from the unstable de 
Sitter node towards the singularity. The region O, whose boundaries are the separatrices coming out from the saddle 
and going to the saddle in the region x > 0, is covered by the trajectories of closed models which begins its expansion 
from the initial singularity, reach a maximal size and then recollapse to the final singularity. The trajectories situated 
over the upper branch of the trajectory of the flat model k = describe the open models expanding towards stable 
de Sitter universe from the initial singularity. They are called oscillating cosmological models. 

The second class of models, described the cosmological models with a squeezing bounce phase in the cosmological 
scenario, is present in Fig. EI Its phase portrait contains a center - non- hyperbolic critical point whose presence makes 
the system structurally unstable. Both phase portraits in Fig. [5] and [3] are topologically non-equivalent in the sense 
of existence of homeomorphism. In Fig. ^ we present a phase portraits of cosmological models with acceleration in 
all time; however there is no matter dominating phase during its evolution. In Fig. [S] it is presented a sub-case which 
starts with accelerating phase. This model is similar to the models considered in loop quantum cosmology which 
started with acceleration at very beginning. 

In the simplest case we obtain a phase portrait equivalent to that presented on Fig. Notice that all trajectories 
are bouncing and trajectories around the center are oscillating. The accelerating region is situated left of the critical 
point. In Fig.[S]the qualitative behaviour of the universes whose late time evolution is dominated by the cosmological 
constant term is presented. In this case we have two disjoint acceleration areas corresponding to that domain of 
configuration space at which the potential is a decreasing function of the scale factor. 

One can imagine different evolutional scenarios in terms of the potential function (see Fig. 1216( 1 . Because of the 
existence of a bouncing phase which always gives rise to the presence of a non-hyperbolic critical point on the phase 
portrait one can conclude 
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forbidden region 



FIG. 6: The model of an accelerating universe given in terms of the potential function and its phase portrait. The domain 
of acceleration is represented by shaded area. The phase portrait as well as the diagram of the corresponding potential for 
bouncing models H 2 /H 2 — f2 m ,os _m — Q n fiX~ n , n > m, m,n = const. All trajectories represent cosmological models in which 
the phase of expansion is proceeded by the phase of contraction. They are closed for oscillating (O) models appearing for 
positive curvature, or open, representing bouncing flat and open cosmologies with a single bounce phase. On the phase portrait 
we find a single critical point of the center type. Because the presence of a critical point of non-hyperbolic type, the system 
is structurally unstable. The motion of the system is determined in the domain of the configuration space V(a) < 0, i.e. for 
a > dmin in the case presented on the figure. Note that for oscillating models there is an infinite number of transitions from 
acceleration to deceleration epoch. For flat and open cosmologies there is only one such transition. 



• the bounce is not a generic property of the evolutional scenario, 

• structural stability prefers the simplest evolutional scenario in which the deceleration epoch is followed by the 
acceleration phase. 

The dynamical systems with the property of such switching rate of expansion, following the single-well potential are 
generic in the class of all dynamical systems on the plane. 

The presented approach to describe dynamics can be extended to the case of cosmological models with scalar 
field. They play an important role in the quintessence conception. To illustrate this let us consider a homogeneous, 
minimally coupled scalar field on the FRW background. The dynamical effects of this scalar field are equivalent to 
the effects of a perfect fluid with energy density and pressure given in the form 

P<t> = ±e</> 2 + [/(0), 
P4, = - U{<j>), 

where e = +1 for standard scalar field and e — — 1 for phantom scalar field, U(<p) is the potential of a scalar field. 

A construction analogous to the presented above is to use the expression for the effective energy density. Let us 
consider for example a universe filled with perfect fluid with pressure p — wp, w — const and a minimally coupled 

2 

scalar field. Then, we can adopt the standard formula V = — Pc 'g a (if we use conformal time dry = dt/a then 
V = - £a f^) and we obtain 

V = -^p W fia^ 3w - ^P<pa 2 . 

After substituting p^ into the above formula and the shifted kinetic term — j^^ 2 m to the kinetic energy of the system 
(remember that the division into the kinetic and potential parts has a purely conventional character), we obtain a 
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2-dimensional Hamiltonian system in the form 

where in above formula we use units in which 8irG = c = 1 (some authors put AttG/3 = c = 1, then V = —p c sa 2 ), <f> 
is a rescaled function (f> — » 0/vo and 

6 

In the special case of the radiation filled universe (w — 1/3) we obtain Hamiltonian system defined on the zero energy 
level 

k 1 

E = -- + -pr.o = const. 
2 b 

Therefore, in the case of potential U (<fi) = ^m 2 4> 2 + \<f) A we obtain the 2-dimensional potential function ip( a i 0) which 
identify quintessence model. The same approach can be adopted to the case of conformally coupled scalar field a well 
as for complex scalar fields. 

It is interesting to investigate what a class of quintessence models which reduces the dynamics of the universe with 
minimally coupled scalar field $ (evolving in the potential E/($)) to the dynamical system of a Newtonian type with 
the potential parameterized by the scale factor V(a). The quintessence model is determined completely once the 
potential function U of the scalar field is given. Thus energy density and pressure is given by 

PeS = -0 5—. 

On the other hand the dynamical system of a Newtonian type is specified completely once the potential V(a) is 
fixed. Therefore if we assume that $ = $(a(t)) then 

p^ = U>'{a) 2 a 2 H 2 {a) + U^{a)) 

p^ = U>'{a) 2 a 2 H 2 {a)-U^{a)) 

and 

<5>'{a) 2 a 2 H 2 {a)-2U{<5>{a)) 



w x (a) 



<f>'(a) 2 a 2 H 2 (a) +2U($(a)) 



where a 2 H 2 (a) = -2V(o). 

Hence the particledike approach can be extended naturally on the class of phcnomcnological quintessence models 
with the parameterized equation of state by redshift (or the scale factor equivalently) |52l l53j . In this approach it is 
assumed the equation of state px = wx{o)px and wx = wxifl) such that 

Px = Px fi a-^ 1+ ^ = pxs^e^ xW = px.o^e' 3 ^ — ^ 
where wx (a) is the mean of the equation of state in the logarithmic scale, i.e. 

/ wx(a)d(lna) 



w x (a) 



J d In a 



The main motivation of the above assumption is the explanation of cosmic coincidence problem: why do all the 
contributions from the vacuum energy density (the cosmological constant) are comparable with energy density of 
matter? To remove the fine-tuning problem it was proposed a simple power law relation wx(a) — woa a , wq, a = const 
(scalling fluid) [54j. If we consider a scalar field parameterized by the scale factor, the class of possible quintessence 
paths is restricted from definition. Note that cosmography measures only average properties of matter expressed in 
kinematics of H(z) and from (1) we obtain [Hlj 



d*\ 2 _ 
~dt) 



px(l + w x ) (11) 



U(^) = ^px(l-wx) (12) 
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where it is assumed that both the scalar field $ and its potential V($>) depends on time through the scale factor, i.e. 
= $(a(i)), [/(<£>(£)) — U(Q(a)). Then equation ifTTl) can be rewritten to the new form 



d§ 
da 



a 2 H 2 



px(a)[l + wx(a)] 



or in terms of redshift z 



d<i> 

dz 



px{z)[l + w x ] 
' H 2 (z)[l + z} 2 



$ = $(z) 



(13) 



(14) 



where we use the standard relation 1 + z = a 1 and ao = 1 is the value of the scale factor at present epoch. The 
potential V in the dynamical system under consideration and the potential of the scalar field can be written as 



V(a(z)) 



a 

~6 



or 



and 



Hence 



V{a{z)) 



1 r 
6 



Pm,oa~ 3 + px,oa-' 3 e- 3 ^ a ™x(a)d(l na ) 
U{a) = p x , a- 3 e- 3 ^ Wx ^ d{Xna Hl - w x ) 



V ( a ) = - g [Pm.oa -3 + (1 - wx^Uia)] a 2 . 



(15) 



(16) 



(17) 



(18) 



Let us return to the general Hamiltonian (|10fl then after substitution (|13|) to the potential V = —p c sa 2 /6, p e g 

P# + Pm 



px(a)(l + w x {a)) 



U(a) 



1 -i 1 2 
6 6 



px{a){l+w x {a)) , 1 n s' 
2 + 2 Pxi - ~ x > 



1 
6 



-\px{a)a 2 - ^Pm.oa 1 = -^Peffa 2 
6 6 6 



(19) 



This means that after parameterization by the scale factor the scalar field quintessence model has the potential 
in the form just prescribed in the particle like description of quintessential models in which the phase space is two- 
dimensional instead of 4-dimensional (<!>, $,a,d) |56l l57|. Let us note that in this case the potential reconstruction 
performed by Rahvar and Mohaved is equivalent to the reconstruction of the corresponding potential function of the 
dynamical system V(a). All dynamics is now squeezed to the plane (a, a). 



III. THE VALUE OF TRANSITION REDSHIFT AND THE HUBBLE PARAMETER AT THE 

TRANSITION EPOCH 

Due to the existence of a simple relation between the luminosity distance d^z) as a function of redshift z and the 
Hubble function 



H(z) = 



d 

dz 



l + z 



(20) 



it is possible to reconstruct the potential function V(a). Note that even if the luminosity distance i obtained accurately 
the potential cannot be determined uniquely as it depends on fi^o - an additional curvature parameter. 

If we assume that the model is flat, then the trajectories of the corresponding Hamiltonian system lies on the zero 
energy level and 



l + z 



(21) 
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FIG. 7: The potential function of the system reconstructed from SNIa data (the Riess et al. sample). 

The result of this reconstruction from the Riess sample is shown in Fig. From this figure we obtain the shape 
of the potential function like the one of the ACDM model. Hence the corresponding phase portrait contains a single 
maximum. The value of redshift at this point we call the transition redshift and denote zt- It is possible to obtain 
not only the qualitative shape of the potential function (an inverted single- well potential) but also some quantitative 
attributes of the model. Especially it is well to know when the switch from the deceleration to the acceleration epoch 
occurs in different cosmological scenarios. The result of fitting of the Newtonian type system with the simple potential 
function given in the simple polynomial form is presented in Table [H] From this estimation we obtain zt — 0.5 which 
can be interpreted as the value of redshift zt at which the switch between deceleration and acceleration epochs takes 
place. 

In Fig.0one can see that the universe is accelerating in the redshift interval in which the potential function decreases 
in respect to z. On the other hand the acceleration takes place in the region where the strong energy condition is 
violated. The idea of testing of energy condition through the measurement of distant SNIa was done by Santos et al. 

M 

Due to the particle like description dynamics of accelerating models one can also estimate the value of the Hubble 
function at the moment of the transition redshift, namely 

H 2 {z T ) = -2(1 + z T fHlV{z T ). (22) 

From equation (|22(l we obtain that 



H(z T ) ~ 0.9H ( 



(23) 
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There arises a basic problem in connection with the above coincidence: Why the value of H(zt) and the present value 
of the Hubble function Hn are of the same order of magnitude? The nature of this problem is similar to the cosmic 
coincidence conundrum (58j |. 

From the definition of the potential function we can obtain a simple interpretation of acceleration of the Universe. 
For example, the deceleration parameter at the present epoch is the slope of the potential (qo = (dV/da) a= i), jerk 
j = "da 1 jh? is related to the convexity of V(x) : jo — — (dV/da 2 ) a =i- Our results also show the advantage of 

using V{x) instead of the coefficient of the equation of state wx{z) to probe the variation of dark energy. 

It have been recently shown the advantages of using px{z) and H(z), instead of the coefficient of equation of state 
[(Olllil- Note that both potential approaches are useful to differentiate between different dark energy propositions. 

The potential function in the neighbourhood of its maximum can be approximated by 



V(a) = V(a T ) + 



d 2 V\ (a-a T ) 2 



da 2 J 2 

Hence we can calculate the energy density near the transition epoch a = ax 

6, 2 , wf d 2 V \ n fd 2 V\ a T fd 2 V 

where (^rjr) = h a Tj^ \a T (p + 3p). The first and second terms are positive and correspond to 1-dimensional 

\ / a—a-j- 

topological defects and positive cosmological constant terms. The third term is negative and scales like 2-dimensional 
topological defects. 

It is useful to define space time metric of the non flat universe in the new variables 



d,s = a 



2 



d(lna) 
7 = 2T 



d^ 



ds 2 =a 2 {dr 2 -d^}, 

where dr = ~P=fpp ■ Therefore, the only nontrivial metric function in a FRW cosmology is the function of V(a) and 
the value of curvature which is encoded in the spatial part of the line element. Hence one can conclude any kind of 
observation based on geometry (cosmography) will allow us to determine a single potential function (for comparison 
see 0). As argued by Padmanabhan, this function is insufficient to describe matter content of the universe and some 
additional input is still required. 

Let us consider now the general properties of dark energy dynamics in terms of the potential. Due to particle-like 
description of cosmology with dark energy, the methods of qualitative analysis of differential equations can be naturally 
adopted. The main advantage of this method is the possibility of investigating all admissible evolutional paths for all 
initial conditions in the geometrical way - on the phase plane. The structure of the phase space is organized by singular 
solutions of the system which are represented by critical points (points of the phase plane for which right hand sides 
of the system vanishes) and phase curves connecting them. Following the Hartman-Grobman theorem the behaviour 
of the trajectories near the critical points is equivalent to the trajectories of linearized system at this point. Therefore 
for constructing of the picture of global dynamic called a phase portrait it is necessary to investigate all critical 
points and their type (determine the stability). The critical points (xq,0) correspond to yo = and (dV/dx) XQ = 0. 
They are saddle point if V xx (xq) < and then eigenvalues of the linearization matrix Aj.,2 = Vxx(xo) are real 

of opposite signs or centres if V xx {xq) > and then eigenvalues are purely imaginary A1.2 = ±iy / V xx (xo)- The 
exceptional case of V xx (xq) = is degenerate. Because only static critical point are admissible due to the constraint 
condition (Friedmann first integral) we obtain that V(xq) = ^k,o- Therefore, because V(x) < 0, the critical points 
are admissible for Sl^.o < 0, i.e. for the closed model only. 

The system linearized around the critical point xq has the form 

(x-x o y=(y-0), (24) 

(y-yoJ= (-^) _ (x-xo), (25) 

which is equivalent to a single differential equation of the second order 
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The solution of the above linear system which approximates the behaviour of the trajectories near the critical point 
is 



where xq = x(to) and xo can be determined from the Friedmann first integral xq = y— 2V{ao) xq . The special choice 
to = ±00 corresponds to the separatrices in and out-going from the static critical point. 

Let us consider now the case of a saddle point. Then the angle of slopes of the separatrices at this point 



tan — 



a fd 2 V 



dx 2 



where a is the angle between eigenvectors at the critical point of the saddle type. 

On the other hand from the reconstructed form of the potential function one can determine the dynamics on the 
phase plane without any information about the value of Qk,o- Hence, we can establish the angle a. It is strictly 
related to the value of jerk because of the relation j — ~(d 2 V/dx 2 ) XQ . Finally we obtain 

a — 2 arctanjo- 

The universe is accelerating in such a domain of the configuration space in which V{x) is a decreasing function of 
its argument. One can calculate the average value of time which trajectories spend during the loitering epoch (a close 
to zero) 

1 f Xn+Ax xdx 1 f X0+Ax dx 

At 



2x Jx -Ax y/-2V(X) 2x Jx -Ax H(x) ' 

where xq is the value of the scale factor at the transition epoch expressed in the units of its present value. 

Hence the time which the model spends in the loitering epoch depends on the transition epoch (z tr ) and the 
preassumed value of Ax - which measures deviation from this stage. For the FRW model with the A term one can 
find exact forms of function At(xo, Ax) in terms of the Jacobi elliptic functions. 



IV. THE GENERIC AND NON-GENERIC GLOBAL EVOLUTIONAL PATHS IN THE MULTI VERSE 

OF ACCELERATING MODELS 

The analysis of full dynamical behaviour of trajectories requires the study of the behaviour of trajectories at infinity. 
It can be performed by means of the Poincare sphere construction. In this approach we project the trajectories from 
center of the unit sphere S 2 — {(X, Y, Z) G R 3 : X 2 + Y 2 + Z 2 = 1} onto the (x, y) plane tangent to S* 2 at either the 
north or south pole (see |50|, p. 265]). Due to this central projection (introduced by Poincare) the critical points at 
infinity are spread out along the equator. Therefore if we project the upper hemisphere S* 2 onto the (x, y) plane of 
dynamical system of the Newtonian type, then 



or 




There is a simple way to introduce the metric in the space of all dynamical systems on the compactified plane. 
If / G C 1 (A4) where M. is an open subset of R n , then the C 1 norm of / can be introduced in a standard way 

||/|| 1 -sup|/(x)| + supp/(x)||, (28) 

where | . . . | and || . . . || denotes the Euclidean norm in R™ and the usual norm of the Jacobi matrix Df(x), respectively. 
It is well known that the set of vectors field bounded in the C 1 norm forms a Banach space (see |5(j, P- 312]). 
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It is natural to use the denned norm to measure the distance between any two dynamical systems of the multiverse. 
If we consider some compact subset JC of A4 then the C 1 norm of vector field / on JC can be defined as 

||/||i = max | /(x) | +max||D/(a;)|! < oo. (29) 

Let E = R™ then the e-perturbation of / is the function g 6 C l {M.) form which ||/ — g\\ < e. 

The introduced language is suitable to reformulate the idea of structural stability given by Andronov and Pontryagin. 
The intuition is that / should be structurally stable vector field if for any vector field g near /, the vector fields / 
and g are topologically equivalent. A vector field / € C 1 (A / I) is said to be structurally stable if there is an e > 
such that for all g E C 1 (A'J) with ||/ — g\\i < e, f and g are topologically equivalent on open subsets of R™ called 
M.. Note that to show that system is not structurally stable on R™ it is sufficient to show that / is not structurally 
stable on some compact JC with nonempty interior. 

It was originally a wide spread opinion that structural stability was a typical attribute of any dynamical system 
modelling adequately a physical situation. The 2-dimensional case is distinguished by the fact that Peixoto's theorem 
gives the complete characterization of the structurally stable systems on any compact 2-dimensional space and asserts 
that they form an open and dense subsets in the space of all dynamical systems on the plane. 

Let us apply this framework to the multiverse of accelerating models represented in terms of a 2-dimensional 
dynamical system of the Newtonian type on the Poincare sphere S 2 . Let the potential function be given in the 
polynomial form. If we assume that the coefficient of the equation of state can be expanded around the present epoch 
(a = 1 or z = 0) in the Taylor series, i.e. 



Px = w x {a)p x , 

w x(a) =E"o w >( 1 - a ) 1 



(30) 



then we obtain from the conservation condition relation 



p X =Px,oa- 3 ( 1+ ^"o(- 1 ) S -)exp{3^(-l) fe ii^^(-irV}. (31) 

k=0 i=k 

Hence if we expand exp{. . . } in formula 1)31(1 then we obtain px as well as V(a) in the polynomial form 



N „ N 



Px = Px,oa-^ + ^o(-^) ) 1 + 3 £(-l) fc <i^L £(-1)^ + I 



k=0 



N „ , k N - r 



i—k n 



k 

k=Q 



(32) 



and 



p m a 2 pxa 2 



V{a) = -^f- (33) 
In the simplest case of w(a) linearized around a — 1 (w = Wo + (1 — a)Wi), we obtain the density of dark energy 
px - Px,oa- 3{1+Wo ~ Wl) exp [3t«i(a - 1)] - P x,oa- 3{1+w °~ Wl) {l + 3w x (a - 1) + 9w 2 (a - l) 2 + • • • } , (34) 

where a -1 = 1 + z, x = a, ao = 1. 

If we consider some subclass of dark energy models described by the vector field [y, — (dV/dx)] T on the Poincare 
sphere, then the right hand sides of the corresponding dynamical systems are of the polynomial form of degree m. 
Then / is structurally stable iff (i) the number of critical points and limit cycles is finite and each critical point is 
hyperbolic - therefore a saddle point in finite domain, (ii) there are no trajectories connecting saddle points. It is 
important that if the polynomial vector field / is structurally stable on the Poincare sphere S 2 then the corresponding 
polynomial vector field [y, — (dV/dx)] T is structurally stable on R 2 f [5(l p. 322]). Following Peixoto's theorem the 
structural stability is a generic property of C 1 vector fields on a compact two-dimensional differentiable manifold Ai . 
If a vector field / g C 1 (A4) is not structurally stable it belongs to the bifurcation set C 1 (A4). For such systems their 
global phase portrait changes as vector field passes through a point in the bifurcation set. 

Therefore, in the class of dynamical systems on the compact manifold, the structurally stable systems are typical 
(generic) whereas structurally unstable are rather exceptional. In science modelling, both types of systems are used. 
While the structurally stable models describe "stable configuration" structurally unstable model can describe fragile 
physical situation which require fine tuning [Isj . 
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In Figs. |HHUEl we show the phase portraits of different evolutional scenarios offered by different propositions of 
solving the cosmological problem. Among different models only the ACDM model (Fig. |SJ (or phantom cosmology) 
and bouncing cosmology in Fig. ITIlk give rise to structurally stable evolutional paths. 

The case of the phantom cosmology requires the additional comment. Let us consider the phantom cosmology with 
wx = 4/3 and dust matter. At the finite domain of the phase plane the system is described by 

x = y (35a) 
1 3 

y = -^m,ox" + 2^ph,o£ 2 (35b) 



and possesses the first integral in the form 



where 



^y 2 + V(x) = -fi fe>0 



V(x) = -^m,ox 1 - ^ P h,oa; 3 - 

To investigate the behaviour of trajectories at infinity (i.e., on the circle at infinity x 2 + y 2 = oo 2 ) it is useful to 
introduce the projective coordinates on the plane. There are two maps which cover the circle at infinity. Let us 
consider one of them 

1 x 

(x,y) h- > (v, w) : v = — ,w = — 

y y 

Then the circle at infinity is covered by y = 0, — oo < w < +oo (for full investigation of the system the second map 
(x, y) i— > (z, u) : z — i, u = - should be studied and both coordinate systems are equivalent if v ^ and 11 / 0). 
In the coordinates (v, w) system (135(1 assumes the form 

^~ = \ (^m,ot' 5 - fl P h,ow i v) (36a) 

dw 1 4 3 c, o , . 

= 2^ni,o« w - -w° + vw (36b) 

where dt/drj — vw 2 . Hence there is one (double) degenerated critical point situated at x = 0, y — oo. 

The naive thinking about the system gives rise to the supposition that the system is structurally unstable because 
of the existence of degenerated critical point in infinity. The objections are the following. All solutions of dynamical 
system x(t,Xo) can be divided on two categories - singular and non-singular. The former are represented by critical 
points in the phase space. The latter are visualized by trajectories joining them. While the trajectories at the finite 
domain of the phase space represent the physical evolution of the evolution, the trajectories at infinity are added to 
the model by the Poincare sphere construction and they cannot represent the physical solution. 

If we go along the trajectories in the future then at some moment of time, say t = tfi na i become to tangent to the 
circle at infinity. This state is called the big-rip singularity. On the circle at infinity we can find a trajectory joining 
the big-rip singularity with the degenerate critical point (x,y) = (0, +oo). Note that this trajectory has no physical 
interpretation. 

In the definition of structural stability itself it appears the assumption that the boundary of the domain at which 
the system is considered is a "cycle without contact" , i.e. there is a simple smooth curve C which does not intersect 
the boundary. Although this assumption bounds a class of systems it makes this sense of notion of structural stability 
to be simpler |63|. 

The parameterization of time, applied for establishing the qualitative equivalence of phantom model with the ACDM 
model, is not of course a diffeomorphism. However, the phantom system in the new parameterization can be treated as 
the model of physical reality. In this case the big-rip singularity is reached for infinity time due to reparametrization of 
time (see Fig.|5J). In this case there is no non-physical trajectory lying on the circle in the infinity. This transformation 
prolongs incomplete trajectories to infinity and the big-rip singularity is a global attractor. From the physical point 
of view this parameterization of the phantom cosmology seems to be more adequate. Therefore, we claim that the 
phantom cosmology is structural stable. However, we must remember that the existence of topological equivalence 
of trajectories of the phantom and ACDM models does not mean their physical equivalence. It is manifested by 
the non-diffeomorphic reparameterization of time. It is an example that the scientific modelling we can use both 
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FIG. 8: The phase portrait of the ACDM model on the projective phase plane RP 2 . Using the Poincare sphere construction 
we represent the dynamics of ACDM models on the compactified plane with the circle at infinity. Note that critical points 
(hyperbolic) at infinity are structurally stable. 

the fragile and structurally stable models. But some of them seems to be more adequate of description of physical 
processes. 

The cosmological models with dynamics as presented in Fig. HOfa , describe for example models in which instead of 
the initial singularity there is a characteristic bouncing phase. There are different situations which realize this typ e of 
evolution: the cosmological models with spinning fluid [fjij or the metric affine gravity (MAG) inspired model |33llr35| . 
All trajectories in Fig. 110b represent bouncing models. 

The qualitative behaviour of extended bouncing models (see Table with the cosmological constant is illustrated 
in Fig. 110b . Because of the presence of an additional center on the phase portrait the models lost the property 
of structural stability which is the attribute of the models in Fig. HOfa . Therefore, they are non-generic in the 
multiverse. An analogous type of evolution appeared as a phenomenological implication of discreteness in Loop 
Quantum Cosmology |6fiq . In general, if we consider a squeezing bounce phase in an evolutional scenario, we obtain 
a fragile model which is structurally unstable. 

Recently Ashtekhar (Loops '05 Conference presentation) suggested that quantum gravity (geometry) can serve as 
a bridge between vast space time regions which are classically unrelated, i.e. in Loop Quantum Gravity, singularity 
is a transitional phenomenon. Therefore, the resolution of the singularity problem of general relativity is replaced by 
approach to singularity with a bounce generated by quantum effects. 

From our point of view such types of evolutional scenarios are not typical in the space of all evolutional paths on 
the plane. 

Moreover, brane world models allow for a transient acceleration of the Universe which is preceded and followed by 
matter domination epoch (deceleration epoch). They admit so-called "quiescent" cosmological singularities [67|], at 
which the density, pressure and Hubble parameter remain finite, while all invariants of the Riemann tensor diverge to 
infinity within a finite interval of cosmic time. From our consideration such models are exceptional, i.e. while time-like 
extra dimension can avoid cosmological singularity by bounce this proposition is not generic in the multiverse. 

The Sobolev metric introduced in the multiverse of dark energy models can be used to measure how far different 
cosmological model with dark energy are to the canonical ACDM model. For this aim let us consider a different dark 
energy models with dust matter and dark energy. We also for simplicity of presentation assume for all models have the 
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FIG. 9: The phase portrait for phantom (wx < — 1) cosmology. The critical point at the circle x = 0, x = oo is degenerate 
but after redefinition of the positional variable x x = x 3 ^ 2 (x \— > x = a; _ ( 1 + 3 "'x)/ 2 m general) and then reparametrization of 
time following the rule r i— * r]: |a; 1 ' 3 dT = dr/, we obtain a "regularized" system for phantom cosmology which is topologically 
equivalent to that in Fig. |H] It is important for our aims that this physically equivalent dynamical system is typical in the 
multiverse. 
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same value of f2 m) o parameters which can be obtained, for example, from independent extragalactic measurements. 
Then, the distance between any two cosmological model, say model "1" and model "2" is 

d(l, 2) = max{\V lx - V 2x \, \V lxx - V 2xx \}, 

where we assumed the same value of the Hq parameter measured at the present epoch for all cosmological models 
which we compare, V\ and V% and their derivatives are only the parts of the potentials without the matter term. 

Because (Vi — Vi) x measures the difference of slopes at any x, the distance d in general becomes a function of x. 
However, for fitting or constraining a model's parameters we use actual observational data obtained at x = 1 epoch. 
Therefore, a closed region C on which we compare predictions of different theoretical models can be chosen such that 
C = {1} x {y : < y < y/—2V(l), V(l) = —1/2}. Of course, it is a closed domain. From the definition of the metric 
d, two models are close if the slopes of tangents at the present epoch are close. This metric can also be expressed in 
dimensionless parameters, then we simply obtain the deceleration parameter instead of slopes. Also, instead of C 1 
metric C r metric can be defined. Then C r metric measures the distance between any two models more precisely by 
comparing their additional acceleration indicators like jerk, snap, crackle or higher derivatives of potential functions. 
In general, the metric is a function of the model's parameters. Therefore with estimating the values of parameters, one 
can immediately determine the distance between any two models - elements of an ensemble of dark energy models. 
Because at present we can only measure the second derivative of scale factor as an acceleration indicator, C 1 metric 
seems to be sufficient to order variety of dark energy models for answering how far the model is from the ACDM 
one and any other. In this way we obtain a ranking of cosmological models from the point of view of closeness to 
the concordance ACDM model. Different models belong to the open ball at the center of which the ACDM model is 
located. One can show the existence of the following inclusion relation [63, l69| 

ci _ ^(lUxfitted) Qtux=— 4/3 q 

"ACDM <- ^phantom e D phantom L J Cardassian ' 1 1 

It is interesting that the above ordering is correlated with one obtained from the Bayesian information criterion 

MM. 



V. CONCLUSIONS 



The main goal of this paper was to investigate the structure of the space of all FRW models which offer the possibility 
of explaining SN la data. To this aim we presented a unified language of dynamical systems of the Newtonian type in 
which the potential function determines all properties of the system. We defined the space of all dynamical systems, 
called the multiverse, of accelerating models. This space can be naturally equipped with the structure of Banach 
space which measures the distance between two models. This metric can be used to measure how far different model 
are from the concordance ACDM model. 

The complexity of different models can be defined in terms of the potential function. This function determines the 
domains of the configuration space in which the Universe accelerates or decelerates. The concept of structural stability 
was used to distinguish generic models in the multiverse. Following Peixoto's theorem we called them typical in the 
multiverse because they form open and dense subsets. The structurally unstable models are exceptional and form 
sets of zero measure in the multiverse. Our main result is that the structural stability property uniquely determines 
the shape of the potential. It is shown that genericity favours inverted single-well shape of the potential function. 

On the other hand, this function could be reconstructed (modulo curvature) from distant type la supernovae data 
and the obtained function is equivalent to the potential function for the ACDM model. The evolutional scenario 
in which the acceleration epoch is proceeded by deceleration is uniquely distinguished. Therefore, while different 
theoretically allowed evolutional scenarios have been proposed and formulated in terms of the potential function, the 
simplicity of the evolutional scenario is the best guide to our Universe - the inverted single- well potential function is 
preferred. In other words, our Universe shows evidence of complexity and at the same time great simplicity which 
allows us to probe its properties with the help of simple models. 

In explanation and understanding of the observational data we use the models (an idea of cosmological models). 
However, between the physical subjects and theoretical notions there is a 1-1 correspondence preserving some relations 
(isomorphism). 

If we require the property of structural stability of the model we implicitly apply what McMullin called cosmogonic 
indifference principle |72| . While the anthropic like principles concentrated on explaining some properties of the 
Universe by the specially chosen model parameter, initial conditions or laws of physics, the indifference type of 
explanation concentrate on searching upon very generic initial conditions and laws of physics which act to produce the 
special configuration (see |jj p. 18]). Stoeger, Ellis and Kirchner argue that indifference principle is more interesting 
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from the physical point of view relative to anthropic type of explanation which is more useful in philosophy. The 
authors claim that the problem of choosing between two principles: special or generic model has a rather philosophical 
(or epistemological) character to which the interpretation of our results is strictly related. 

Therefore if we concentrate on searching for a very generic class of dark energy models or modification of the 
FRW equation that produce the special configuration we now enjoy - an accelerating Universe, then the postulate of 
structural stability naturally gives rise to such a situation. 

A natural question is whether the inverted single-well potential is favoured over other more complex models with 
double, triple etc. accelerating phases. This can be addressed by using the Bayesian information criteria (B1C) 
of model selection. Our analysis confirms that there is no strong reason for inclusion of extra complexity (more 
accelerating epochs) and the model with a single acceleration epoch is favoured over the others by supernovae data 

Different conclusions can be made from our analysis but the answer to the question put in the title seems to be 
especially tempting. We are living in an accelerating Universe because simplicity is the best guide to our Universe. 
Moreover, because our Universe is typical (generic) and its model can be discovered by using the approximation 
method starting from a simple (possibly naive) model - the ACDM (or phantoms). 
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